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ABSTRACT 

We report on an 105 ks Suzaku observation of the supernova remnant CTA 1 
(G119.5+10.2). The Suzaku soft X-ray observation was carried out with both tim- 
ing mode and imaging mode. A «~ 10' extended feature, which is interpreted as a 
bow-shock component of the pulsar wind nebula (PWN), is revealed in this deep ob- 
servation for the first time. The nebular spectrum can be modelled by a power-law 
with a photon index of ~ 1.8 which suggests a slow synchrotron cooling scenario. The 
photon index is approximately constant across this extended feature. We compare and 
discuss our observations of this complex nebula with previous X-ray investigations. We 
do not obtain any significant pulsation from the central pulsar in the soft (0.2-12 keV) 
and hard (10-60 keV) X-ray data. The non-detection is mainly due to the loss of the 
precise imaging ability to accurately determine the source contribution. The spectra of 
XIS and HXD can be directly connected without a significant spectral break according 
to our analysis. Future observations of NuSTAR and A stro-H would be able to resolve 
the contamination and provide an accurate hard X-ray measurement of CTA 1. 
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1 INTRODUCTION 

CTA 1 (G119.5+10.2) is a composite super nova remnant 
(SNR) with centrally b right X-ray emission (| Seward et al.l 
Il995l : lsiane et alj|l997h . This SNR has a radio shell with a 
radius of ~ 54', and its kinematic distance estimated from 
the a ssociated Hi emission is 1.4 ± 0.3 kpc |Pineault et al.l 
Il993h . The X-rays fro m CTA 1 comprise both thermal and 
non-thermal emission jSeward et al.|[l995l h While the outer 
regions mainly consist of shock-heated thermal emission, 
the non-thermal component dominates the central emission, 
which has led to the sp eculation of a hidden pulsar in this 
SNR (|Slaneet all [199% . In addition to the radio and X- 
ray observations, the launch of Compton Gamma-Ray Ob- 
servatory (CGRO) provided a further insight on the na- 
ture of CTA 1. Energetic Gamma Ray Experiment Telescope 
(EGRET) on board CGRO has reveal ed a 7-ray source, 2EG 
J0008+7307 ^Thompson et alj Il995l ). in the 2nd EGRET 
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catalo gue (or 3EG J0010+ 7309 in the 3th EGRET cata- 
logue; |Hartma^et^l][l99^), which was suggested to be as- 
sociated with CTA 1 in view of the positional coincidence 
|Brazier et al.|[l99^ ) . The lack of variability of this extended 
7-ray source as well as its spectral resemblance to a typical 
7-ray pulsar further suggested that a neutron star should re- 
side in this SNR. Nevertheless, the limited photon statistics 
obtained by EGRET did not allow any meaningful pulsation 
search in the MeV— GeV regime. 

Observations with ROSAT have revealed an X-ray 
oint source, RX J0007.0+7302, at the center of CTA 1 



Seward et all Il995h . It was also identified as the X- 
ray counterpart of 2EG J0008+7307 / 3EG J0010+7309 
jBrazier et all 1 1998T ) . This point source was identified as a 
neutron star candidate with the high 7-ray-to- X-ray and 
X-ray-to- optical flux ratios as well as its X-ray spectral be- 
haviours (|Slane et al.lll997h . which were found to be similar 
to those of the only one radio-quiet 7 - ray pulsar known at 
that time — Geminga (|Brazier et al 1 1 19981 ; lHalpern et al.l 
2004). Utilizing the high resolution imaging capability of 
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Chandra, lHalpern et alj (|2004l ) have uncovered a torus+jet 
feature, which resembles the morphology of a pulsar wind 
nebula (PWN) as seen in other bright 7-ray pulsars. Further- 
more, the sub-arcsecond angular resolution of the Chandra 
image provided a precise position, which can facilitate multi- 
wavelength investigations. To identify the neutron star na- 
ture unambiguously, several studies att empted to search for 



the spin period of RX J0007.0+7302 dBiggs fc Lvn 



Nice fe Saver Il997l ; ISlane et alJ 120041 : lHalpern et al 



1096 



2004 



Lin fe Chandl2005f ) in different energy bands. Nevertheless, 



no conclusive evidence for pulsations in radio/X-ray was 
yielded in these previous studies. 

With the unprecedented sensitivity of the Large Area 
Telescope (LAT) onboard Fermi Gamma-ray Space Tele- 
scope, a 7-ray pulsar, LA T PSR J0007+73 03, was even- 
tually detected in CTA 1 (jAbdo et al]|2008l ). This source 
has a spin-period of ~ 315.86 ms and was the first radio- 
quiet 7-ray pulsar detected by LAT. Its spin parameters 
imply a characteristic age of ~ 14000 yr, which is consis- 
tent with the kin ematic age of CTA 1 (~ 10000 - 20000 yr; 
ISlane et al.lll997h . This strongly suggests an intrinsic as- 
sociation between the pulsar and SNR. With the discov- 
eries of more radio-quiet 7-ray pulsars |Abdo et al.l 12009] ; 
ISaz Parkinson et al.l[20ich . we can now examine them as a 
unique class. It has been speculated that the main distinc- 
tion between radio-quiet and radio-loud 7-ray pulsars is sim- 
ply a geometrical effect. In order to further investigate their 
emitting properties, X-ray observations are deeply required. 

Although the periodic signals of these 7-ray pulsars can 
be significantly detected, Geminga is the only one radio- 
quiet 7-ray pulsar with a detected X-ray pulsation before 
2009 (|Halpern fc Hordll992f ). In order to investigate the pos- 
sible origin of the pulsed X-rays from the radio-quiet 7-ray 
pulsars as a whole class, searches for the possible pulsations 
from these pulsars have been conveyed with dedicated X-ray 
observations, including PSR J0007+7303. With a ~ 120 ks 
observation with XMM-Newton (hereafter XMM) , the X-ray 
pulsation of PSR J0007+ 7303 have been independe ntly de- 
tected bv lLin et all l|201Ch and lCaraveo et alJ lj20ld ). Apart 
from the X-ray pulsation, this deep XMM observation also 
reveals the nebular emission extending for ~ 2' in south-east 
direc tion away from the pulsar (see Fig. 1 in lCaraveo et al.l 
20101'). Togeth e r wit h the torus and jet feature found by 



Halpern et all | |2004 ). all previous X-ray investigations sug- 



gest LAT PSR J0007+7303 has a complex PWN structure. 
To continue the investigation, we have observed this inter- 
esting object with Suzaku, which enables a detailed analysis 
in both soft and hard X-ray bands. In this paper, we report 
the results obtained from analysing the Suzaku data. 



2 OBSERVATIONS 



Suzaku (|Mitsuda et al.l [20071 ) was built with two main pay- 
loads, five ne sted conical thin-foil g razing incidence tele- 
scopes (XRT; [Serlemitsos et al.ll2007t ) with the effective en- 
ergy range of 0.2-12 k eV and a hard X-ray detector (HXD; 
iTakahashi et alj|2007f ) with the effective energy range of 10- 
600 keV. There are two major detectors for spectroscopy 
(X-ray spectrometer and X-ray imaging spectrometer) in 
Suzaku and each of them is set at the focus of an XRT. 
Our observational results on the temporal and spectral anal- 



yses mainly come from the X-ray imaging spectrometer 
fXIS; lKovama et alj|2007h and the HXD. The XIS consists 
of three front-illuminated (FI: XIS0, XIS2, XIS3) and one 
back-illuminated (BI: XIS1) CCD detectors, but only the FI 
chips provide the observations in timing/P-SUM mode. Be- 
cause the time resolution of HXD is 61 fj,s and that of XIS in 
the timing mode is ~ 7.8 ms, our data is precise enough for 
the periodicity search to a radio-quiet 7-ray pulsar with the 
spin period in tens or hundreds of milliseconds. However, 
two of the XIS units with FI chips, XIS0 and XIS2 suffered 
a catastrophic damage on 2006 November and 2009 June re- 
spectively. Since no useful data can be obtained with XIS2 
and a part of the segment in XIS0 is lost, only the observa- 
tion of XIS3 unit can be used for temporal analysis. These 
events seriously decrease the probability to obtain the peri- 
odic signals from the faint sources in the soft X-ray band. 
There are also two independent detector systems used for the 
observations of HXD. The positive-intrinsic-negative (PIN) 
silicon diodes are sensitive below ~ 60 keV in the FOV of 
34' x 34', and the Gadolinium Silicate/Bismuth Germanate 
crystals (GSO/BGO) phoswich counters can detect photons 
above ~ 30 keV in the FOV of 4° .5 x 4° .5. Because the data 
obtained from HXD have no imaging capability, we only did 
the analysis of the HXD-PIN data to avoid the contamina- 
tion caused by the large background size of the HXD-GSO, 
which is much larger than the size of CTA 1. 

Our observation for CTA 1 was carried out on 2010 Jan- 
uary 8 with total exposures of ~ 105 ks in the soft X-ray 
band and ~ 59 ks in the hard X-ray band, and it was aimed 
at the HXD nominal position to RX J0007.0 +7302 with 
(J200 0) R.A -00' l 07 m 02 s .2, decl.=+73°03'07" (jSlane et all 
l2004h . The XIS0 and XIS1 data of Suzaku were observed 
in the normal mode without window option, while the 
Suzaku/XIS3 data observed within 128 rows were com- 
pressed into one dimension to attain a higher time res- 
olution. All the data reduction and the spectral analy- 
ses were performed using XSELECT (ver. 2.4a), FTOOLS 
(ver. 6.9) and XSPEC (ver. 12.6.0) of HEASOFT (ver. 6.9) 
with the latest Suzaku/XIS calibration files (20101108) and 
Suzaku/KXD calibration files (20101202). 

2.1 XISO & XIS1 

The XISO and XIS1 observations of 105 ks exposures were 
divided into two editing modes of 3 x 3 and 5x5 with ~ 85 
and ~ 20 ks, respectively. Because each pixel on the CCD 
was read out every 8 s in the normal mode, the time reso- 
lution of the XISO and XIS1 data is not accurate enough to 
search the periodic signal in the range less than one second. 
However, the XISO and XIS1 data are still adequate for the 
spectral analyses of the X-ray point source and the PWN to 
compare with the results yielded by other X-ray missions. 
The adoption for the center of the point source is referred 
to the aim point of the Suzaku observation. Because RX 
J0007. 0+7302 is seriously contaminated by the surround- 
ing PWN, we extracted the spectra of RX J0007.0+7302 
only within l'.5 central circular region instead of the regular 
choice of 3' circular region for a point source. This adoption 
still contains ~ 75% of the encircled energy fraction (EEF). 
We extracted the background of the pulsar from a concentric 
annulus of radii 100" < r ^ 180" centered at the aim point 
of this Suzaku observation. The spectra extracted from the 
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Figure 1. Vignctting-corrected images of the field around RX J0007. 0+7302 as observed by Suzaku/XISO (panel a) and Suzaku/XISl 
(panel b). A PWN component with an extent of ~ 10' is revealed. The lost of the row image due to a putative micrometcorite hit can 
also clearly be seen in Suzaku/XISO. The inset shows the close -up view of the lxl region (illustrated by the black square) around the 
pulsar as seen by Chandra (adopted from lHalpern et alj|2003 l. For all images, top is north and left is east. 




point source comprise 695 counts for XISO and 401 counts 
for XIS1 after subtracting the background. We also consid- 
ered the point source selection in two other different sizes 
with 2' and 3' in radii using the same center as shown in 
Table [T] Such selections of the point source contain ~ 85% 
and ~ 95% of the EEF, respectively. The background of 
these selected regions were determined within a concentric 
annulus of radii 2' < r < 4', and 3' < r < 4'. 5. Different 
choices for a point source region are used to verify that the 
significant excess contributed by the surrounding PWN can 
be detected with the increasing size. 

For the PWN, we firstly examined its source extent by 
applying adaptive smoothing to both XISO and XIS1 im- 
ages. The images shown in Fig. [T] were generated with tak- 
ing into account the vignetting of the XRT. We used the 
task "xissim" to simulate a flat field image for our purpose, 
and we displayed the obtained image after dividing by the 
flat field. A large structure with an extent of ~ 10' appears 
to extend eastward from the pulsar. In view of the possi- 
ble spectral steepening due to synchrotron cooling, we per- 
formed a spatially-resolved spectral analysis of the PWN. 
The extraction regions are illustrated in Fig. [2] We split the 
PWN into three regions, namely the inner part (a circular 
region within a radius of 2'), the middle part (an annulus 
with radii 2' < r < 4') and the outer part (a panda re- 
gion with radii 4' < r < 8'). The background spectra are 
sampled from the source free low-count regions from each 
detector respectively. We only consider the spectra obtained 
with the 3x3 mode to perform the spectral analysis of the 
PWN because of its better photon statistics. To yield better 
photon statistics to describe the spectral behaviour for the 
PWN, we have rebinned all the spectra before the analysis. 



Figure 2. Illustration of the regions used to extract the 
spectra from different parts of the PWN associated with 
RX J0007.0+7302. 



We generated the spectrum with > 200 cts/bin for the inner 
part, and with > 500 cts/bin for the middle and outer parts. 

We generated the associated response matrix (rmf) and 
auxiliary response (arf) files by the HEASOFT command 
of 'XISRMFGEN' and 'XISSIMARFGEN'. The spectrum 
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obtained from each detector was rebinned with a minimum 
of 30 counts per channel to ensure \ 2 statistics. We then 
introduced a constant in a simultaneous fit to account for 
the cross-calibration mismatch between different detectors. 
We also considered different circular sizes to describe the 
contribution of the point source. Details about the spectral 
results are presented in § 3. 



2.2 XIS3 

The XIS3 data was investigated by the timing/P-SUM mode 
with a long duration of ~ 212 ks and its image was com- 
pressed into one dimension without the information along 
Y-axis to obtain higher time resolution. Our procedures to 
resolve the point source in a one-dimensional image follow 
the stand ard recipe for reduc ing XIS data with the tim- 
ing mode ()Matsuta et alJfeoid ). We also redid the event re- 
processing to consider the grade filtering and generated the 
new cleaned event file using the script 'XISREPRO.XCO' 
with the proper selection criteria. For the P-SUM mode, 
we need to remove the hot pixels on the image manually. 
To identify those hot pixels, we accumulated all the pho- 
ton events along the coordinate of ACTX. We first removed 
those columns that have photons significantly more than 
others, and then we counted the distribution of events in 
the remaining columns and removed those pixels with the 
photons exceeding 3<r than the median of the distribution. 

After excluding all the hot pixels, the source is still too 
marginal to be resolved in the histogram of ACTX versus 
photon events of the XIS3 data. We then tried to rebin 
the distribution of counts along ACTX with several pix- 
els as shown in Fig. [3] A point-like source with 112 pix- 
els (~ 118") can be roughly detected close to the center 
when we rebinned with 32 pixels. There are 58255 pho- 
tons for this source within 0.2-12 keV and the significance 
of this detection is ~ 3.3a. We also checked the absolute 
arrival time of the photon events and determined the time 
required for readout by the ACTX position of the source 
in the XIS1 image. In order to perform a precise tem- 
poral analysis, we corrected the arrival time with the de- 
crease of 31.2 ms for all the source photons. The Solar sys- 
tem barycentric time correction was then implemented with 
the task 'AEBARY CEN' at (J2000) R.A -Q0 ; '07 m 01 s .56, 
decl.=+73°03'08".l (|Halpern et all |2004 ) to yield a time 
list, which can be used to examine the pulsation of 
RX J0007.0+7302/PSR J0007+7303 in the soft X-ray band. 



2.3 HXD 

The HXD observation was conducted in nominal pointing. 
Since the HXD is a non-imaging instrument, we can not 
separate the contribution of the pulsar from the surrounding 
nebula. However, a spectrum can be acquired to describe the 
behaviour of the total hard X-rays emitted from the entire 
CTA 1. 

We only focus on the HXD-PIN detections of the 
hard X-rays emitted from/around CTA 1. In order to 
examine the periodic signal of the pulsar, we also ap- 
plied the Solar system barycentric time correction with 
the task 'AEBARY CEN' at (J2000) R A -Q0 fe 07 m 01 s .56, 
decl.=+73°03'08".l jHalpern et alj|2004 ) on the HXD-PIN 
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Figure 3. Histogram of ACTX versus photons distribution of 
the XIS3 data. The black solid line marks the average photons of 
each pixel after rebinning the ACTX with 16 pixels, and the red 
one shows the distribution with the bin of 32 pixels. The green 
solid line is set as 536 to demonstrate the median of the photons 
distribution. A point-like source can marginally be resolved within 
the range of 496th-608th pixel in the red line. 



observation. The events in the observation were restricted 
according to the effective energy range (10 - 60 keV) of the 
PIN detector, and 26661 photons were obtained for temporal 
analysis after the screening of the energy. 

We also subtracted the non-X-ray background (NXB) 
and cosmic X-ray background (CXB) to obtain the HXD- 
PIN spectrum. The simulated NXB was directly obtained 
from the Suzaku Data Center, and the CXB flux t hat is 
~ 5% of the NXB for PIN could be derived from iBoldtl 
l|l987h . 



CXB(E) = 9.412 x 10" 3 x (E/lkeV)' 1 -'' 

photons cm" 



x exp(-E/40 keV) 



(1) 



These background spectra correspond to epoch "6" for the 
flat PIN response matrix. We also evaluated the contri- 
bution of the Galactic Ridge X-ray Emission (GRXE) us- 
ing the Suzaku observation of a nearby blank sky (OB- 
SID=504039010) at (1, b)=(123°.9, 10°. 0), which is 4.3 de- 
gree apa rt from the CTA 1. Following the same method as 
stated in lEnoto et all (|2010h we assume d a power-law spec- 
trum with a fixed photon index at 2.1 (|Valinia fc Marshall! 
ll99Sf ). The contribution of the GRXE within 14 - 35 keV 
was evaluated to be < 4.4 x 10 -4 photons sec -1 , which only 
corresponds to <0.2% of the NXB. Thus, the contribution 
of GRXE is negligible when we analyse the HXD-PIN spec- 
trum of CTA 1. After the background subtraction, we re- 
binned the spectrum to ensure that the photon numbers in 
each channel are larger than 30. 



3 RESULTS 

3.1 Timing Analysis 

The expected spin frequency of the pulsar in our Suzaku data 
can be inferred fr om the contemporaneous Fermi ephemeris 
(|Rav et aLlkoill ). The start of GTI for our Suzaku/XlS3 
and HXD-PIN data is at epoch MJD 55204.6143077. The 
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Table 1. Best parameters for the spectral fits to the pulsar in 
CTA 1 using Suzaku/XISO and XIS1 observations. 



3.1652 3.1654 3.1656 3.1658 3.166 3.1662 3.1664 




3.1654 3.1655 3.1656 3.1657 3.1658 3.1659 3.166 3.1661 
Trial Frequency (s ) 

Figure 4. Periodicity search of the Suzaku data by the H test. 
Each panel presents the search with 200 independent trials close 
to the expected spin frequency of PSR J0007+7303 at the obser- 
vational epoch. The upper panel shows the searching results of 
the Suzaku/XIS3 data, and the width of each independent trial 
is ~ 4.7 X 10 -6 s . The lower one shows the searching results of 
the Suzafcu/HXD-PIN data, and the width of each independent 
trial is ~ 6.0 X 10" 6 s -1 . 



spin frequency of PSR J0007+7303 corresponding to this 
epoch can be inferred as 3.165750271(6) s" 1 with the first 
frequency derivative of -3.6136(2) x 10" 12 s~ 2 (the effect 
of the second frequency derivative is negligible). We search 
for periodic signals with the known first frequency deriva- 
tive around the predicted freq uency in SuzakufKIS S and 
HXD-PIN data using fl-test (|de Jager et alj 1 19891 ') . The 
Jz~-value is only 11.7 corresponding to the expected spin 
frequency of 3.165750271(6) s" 1 for the Suzaku/XIS3 ob- 
servation and only 6.7 for the Suzaku/HXD-PIN obser- 
vation. Fig. [4] demonstrates the periodicity search of 200 
independent trials around the expected spin frequency of 
PSR J0007+7303. We conclude no significant pulsed detec- 
tion can be yielded from our Suzaku data both in the soft 
and hard X-ray bands. 

Comparing with the soft X-ray pulsation detected 
by th e XMM observation (|Lin et all |2010| ; ICaraveo et all 
2010), our Suzaku data lose 1-dimensional imaging ability 
in the XIS3 timing mode and have no image ability in the 
Suzaku/WXT) investigation. In the Suzaku/XlS3, we can only 
have a marginal detection to resolve the point source with 
~ 58000 photons as shown in section §2.2. However, the 
counts gathered in the source region are much more than 
we expected base d on the spectral behaviour provided by 
ISlane et al.l (|l997h . We noted that most of the additional 
counts obtained from the P-SUM mode data are caused by 
the instrumental background . For the data detected in the 
hard X-ray band, the main contribution of the hard X-ray 
photons is ascribed to the X-ray back g round . If we assume 
the pulsed spectrum of ICaraveo et all (|2010h . the expected 
S/N ratio of the Suzaku investigation is too low (< 2a) to 
yield the periodic signal from the pulsar. 

We also examined this spin frequency using ~ 3-month 
Fermi archive (2009 November 22 - 2010 February 21) on the 
epochs close to our Suzaku observation. The effective 7-ray 
photons of PSR J0007+7303 were restricted in a 1-degree 
circle centered at (J2000) R.A.=1°.7565, deel .=+73°. 05225 





RX 


J0007.0+7302/PL 


Parameter 


Small 6 


Medium'' 


Large 6 



Nh (cm 2 ; fixed) 


2.8 x 10 21 


2.8 x 10 21 


2.8 x 10 21 


r 


1 76+ ' 27 

L - ,D -0.25 


1.80+0.15 


1 «^+ - 14 
l.»b_ 12 


*JC,PL X lO" 13 


3.7+1.1 


4.6 + 1.0 


7.1 ± 1.3 


(ergs cm -2 s _1 ) a 








X 2 /dof 


39.0/46 


78.2/85 


157.4/143 



a The unabsorbed flux is measured in the energy range of 0.3-10 
keV. 

b 'Small, medium and large' indicates the source selection with 
different size. The small, medium and large point source regions 
were adopted with the same center and circles with 1'.5, 2' and 
3' in radii, which correspond to the EEF of ~ 75%, 85% and 
95%. The background contributions were determined within a 
concentric annulus of radii 100" < r ^ 180", 120" < r ^ 240", 
and 180" < r < 270", respectively. 



(jHalpern et alj l2004h within the energy range from 100 
MeV to 300 GeV. We determined the spin frequency of the 
pulsar with the ff-test using the task of 'GTPSEARCtF 
from the Fermi Science Tools ver. v9rl5p2. At epoch MJD 
55204.6143077, the maximum statistic was obtained at 
3.16575026(1) s _1 with the first frequency derivative of 
—3.6136(2) x 10~ 12 s~ 2 , and this result is consistent with 
the aforementioned prediction as well. 



3.2 Spectral Analysis 

We derived the soft X-ray spectrum from XIS0 and XIS1 
onboard Suzaku. Based on our selected region for the X-ray 
point source, the ratio of net source counts for the pulsar to 
the total photons in the source region considered in the spec- 
trum of XIS0 is more than 25% but in that of XIS1 is only 
~ 10%. This is caused by the fact that the BI chip (XIS1) has 
a higher effective area at low energies. We only considered 
the spectral fits in the range 0.5-10 keV and ignored those 
photons outside this energy range to avoid large uncertain- 
ties. A cross-calibration term was included to correct the dif- 
ference among the XIS spectra, and the result is consistent 
with the fit to the spectrum of any individual detector. The 
absorption was fixe d at 2.8 x 10 21 cm" 2 according t o previ- 
ous measurements (|Slane et al.ll 19971 ; iLin et al]|2010l) , whic h 
is consistent with optical extinction (jHalpern et al]|2004l ). 
A single power-law provides an acceptable fit as shown in 
Table [T] and Fig. [5] and our data are not good enough to 
give a well constraint for an additional model component, 
such as a thermal black body model or a m agnetized neutron 
star atmosphere fNSA; [Zavlin et aLlll996l ) model. However, 
we note that the source fluxes derived from different sizes 
seem larger than those of previous studies and increase with 
source sizes; these excesses might be contributed from the 
PWN. 

On the other hand, our deep Suzaku observation pro- 
vides a desirable opportunity to investigate the PWN, 
particularly for the faint emission with the extent up to 
10' (cf. Fig. [TJ. Previous observations with Chandra 



(|Halpern et al.l(2004T l and XMM (|Caraveo et alj|20ich failed 
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Energy (keV) 



Figure 5. Spectral fits to the power-law model. The first row presents the spectral fits to the point source centered in CTA 1 with a 
selected circle of 1'.5 and 2' in radii, respectively. The second row presents the spectral fit to the point source centered in CTA 1 with 
a selected circle of 3' . The black and green data sets refer to the XISO observations of 3 X 3 and 5x5 modes respectively, while the red 
and blue data sets refer to the XISi observations of 3 X 3 and 5x5 modes. The detailed parameters of these fits are shown in Table [T] 



Table 2. Best-fit parameters of the power-law spectral com- 
ponents inferred from different regions of the PWN around 
PSR J0007+7303. a 



Region 


r 


/ P wn (0.3-10 keV) 
erg cm -2 s -1 


X 2 /d.o.f. 


Inner part 
Middle part 
Outer part 


1.89 ±0.08 
1.79 ±0.09 
1.92 ±0.11 


1.2l° 2 X 10~ 12 
(1.4 ±0.2) X lO" 11 
1.5+°i x 10- 11 


30.0/28 
38.1/34 
38.4/41 



a Column density is fixed at Njj = 2.8 X 10 21 cm~ 2 for all 
these spectral fit. An additional power-law component to ac- 
count for th e pulsar emissi o n wit h the parameters fixed at those 
inferred by ICaraveo et al.l 1 I2OIOI ) (i.e. T = 1.36 with /psr = 
6.5 X 10 -14 erg cm -2 s _1 ) is included with a multiplicative factor 
for scaling the pulsar contributions in each region accordingly (i.e. 
85%, 10% and 5% for inner, middle and outer parts respectively. 

in detecting this large but faint diffuse structure; this is likely 
caused by the relatively small collecting area of Chandra and 
the high instrumental background of XMM respectively. 

In order to properly constrain the non-thermal X-ray 
contributions from the PWN, we also fixed the column den- 
sity at Nh = 2.8 x 10 21 cm -2 . Besides a power- law model 



for PWN emission, we also included an additional power- 
law component to account for the X-rays fr om the pulsar 
with t he parameters fixed at those inferred bv lCaraveo et all 
(|2010t ). For each of the three regions in our consideration, 
we scaled the spectral component of the pulsar with the 
EEF or the contribution from the wing of the PSF cen- 
tered at the pulsar accordingly by multiplying a constant. 
The best-fit parameters of the power-law component of the 
PWN inferred from the spatially-resolved spectral analysis 
are summarized in Table [2] Taking the statistical uncertain- 
ties into consideration, no evidence for spectral steepening 
can be found from this observation. The whole diffuse X-ray 
feature can be modelled with a power-law with a photon 
index of T ~ 1.8. 

We also derived the hard X-ray spectrum from the 
HXD-PIN. After subtracting the background, the net counts 
of the hard X-rays collected for the spectrum is only 6.6% in 
total. Because there are abnormal accumulations for photons 
at both limits of the effective energy boundary in our HXD- 
PIN observation, we only consider the hard X-ray spectrum 
within the energy range of 14 — 30 keV. The spectrum fit- 
ted to a power-law model with the photon index (P) of 3.0 
is shown in Fig. [S] We also checked the systematic uncer- 
tainty of the background from other observations, and the 
NXB accuracy should be within 5% level in the 15-40 keV 
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Figure 6. Spectral fit to the Suzafcu/HXD-PIN observation in 
14 — 30 keV. The solid line represents the best-fit to a single 
power-law with the photon index of 3.0, and corresponds to the 
unabsorbed flux of ~ 8.0 X 10 — 12 ergs cm -2 s _1 in 10-50 keV. 
The red line demonstrates the systematic errors which correspond 
to ~ 5% of the NXB distribution. The bottom panel shows the 
residuals in terms of a. 

around the CTA 1 observation. The nominal uncertainty of 
the adopted background model is also indicated in Fig. [6] 
Because the effective data points of our hard X-ray spec- 
trum are few, the uncertainty is very large and the best 
fit can not be well determined. We fixed the absorption as 
2. 8xl0~ 21 cm" 2 acc ording to the column density obtained 
in lSlane etafl (|l997l ), and a photon index of Y = 3.33±fH 
(90% confidence level) with xl = 1-15 for 10 d.o.f. can be 
yielded in the investigation of our hard X-ray spectrum. 
Taking into account the wide spread of the acceptable pho- 
ton index, the total unabsorbed flux in 10-50 keV for the 
whole FOV of the HXD-PIN inferred from these accept- 
able fits distributes from ~ 7.8 x 10 -12 ergs cm -2 s _1 to 
~ 1.1 x 10 -11 ergs cm -2 s , and this signal stays between ~ 
3-5% level of the NXB in the 90% confidence level. The hard 
X-rays detected in the region around CTA 1 is also consis- 
tent with previous INTEGRAL investigation (|Sturner et al.l 
|2004 ). 



4 DISCUSSION 

We have investigated the central region of CTA 1 with a deep 
Suzaku observation. For the X-ray flux of RX J0007.0+7302 
as estimated by the XIS data, we no te th at it is larger 
than that reported by iLin et ail (|2010t ) and ICaraveo et al.l 
|2010l ) through the XMM observation. This discrepancy can 
be ascribed to the fact that the PSF of Suzaku is consid- 
erably wider than that of XMM. In view of this, the con- 
tribution of the point source is likely to be contaminated 
by the surrounding PWN. It is the main reason that we 
can not obtain the real spectrum of the pulsar. Because 
of the loss of the imaging ability in the Suzaku data for 
the periodicity examination, we can not separate the instru- 
mental effect in the soft X-ray band and the serious con- 
tamination from the X-ray background in the hard X-ray 



band. These problems prevent the pulsed detection from 
the current Suzaku observation. Since the Suzaku data do 
not provide a constrained result for the timing properties of 
RX J0007.0+7302/PSR J0007+7303, we can not obtain the 
phase-resolved spectroscopy and will not discuss the nature 
of this point source any further in the following statements. 

We have also examined the properties of the PWN with 
the XIS data. A single power-law can provide a good fit and 
we did not detect any additional thermal X-rays emitted in 
the central region of CTA 1 which is co nsistent with previ- 
ous investigation (e.g. ISlane et al.lll997l ). Both the flux and 
the photon index obtained in our independent investigation 
are fully co nsistent with that in ferred by the recent XMM 
observation (|Caraveo et alj|20ich . 

We would like to discuss the possible nature of the neb- 
ula that extends eastward as discovered by our deep Suzaku 
exposure. The feature revealed by Suzaku extends as large as 
~ 10' while the XMM observations only detected an nebu- 
lar feature with an extent of a few arc-minute, which might 
correspond to th e relative brighter co mponent of the neb- 
ula (cf. Fig. 1 in[c aravco et alj|2010h . The spectral steep- 
ness and the flux of the inner part in our investigation are 
consiste nt with their results (refer to "Outer PWN" in Ta- 
ble 1 of ICaraveo et alJl201Ch . The morphology of this ex- 
tended feature appears to be asymmetric (see Fig. [TJ which 
resembles those of bow -shock nebulae, e.g. PSR J1747-2958 
|Gaensler et al.l 120041 ') . Such extended feature are usually 
along the direction of pulsar motion and behind the bow- 
shock, therefore they are interpreted as the synchrotron ra- 
diation from the flow of particles coming out of the PWNe. 

It is instructive to discuss the plausible relation between 
this large feature and the torus+jet of a smaller scale re- 
ported in Figure 1 of Halpern et al. (2004; also shown in 
the inset of our Fig. [TJ. The Chandra image shows that the 
jet have an extent of ~ 16 toward south and bends to 
southwest at the far end from the pulsar (see Fig. [TJ. On 
the other hand, the compact torus-like feature with a ra- 
dius of ~ 3 around the pulsar appe ars to be elongated in 
a direction perpendicular with the iet. lHalpern et all (|2004l ) 
interpreted this system as an equatorial torus with the jet 
emitted along the rotation axis, which is observed in many 
PWN systems (cf. |Ng fc Romani|[20o3 for a review). 

If the ~ 10 feature is indeed resulted from the bow- 
shock as aforementioned, this suggests the pulsar is possibly 
moving towards west. This speculated pulsar motion is al- 
most perpendicular to the portion of the jet close to the 
pulsar (see Fig. [JJ. The origin of the velocity of pulsars is 
still an open question. One possible mechanism is the asym- 
metric supernova explosion. Apart from giving rise to the 
kick velocity, it also contributes to the initia l spin of a pul- 
sar (|Spruit fc Phinnevl[l998l ; lLai et al] 120011 ). If the initial 
spin angular momentum of RX J0007. 0+7302 is dominated 
by the kick process, the jet emitted along the rotational 
axis should be more or less perpendicular to the direction 
of proper motion, which can explain the relative orientation 
between jet and the ~ 10' feature. Therefore, measuring the 
proper motio n of RX J0007. 0+730 2 with multi-epoch X-ray 
imaging (e.g. iHui fc Becker! 120061 ) and/or dedicated 7— ray 
pulsar timing solution yielded from LAT data can help to 
constrain the relation among various extended feature of this 
complex system. 
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Figure 7. Broadband absorbed spectrum of PSR J0007+7303, 
PWN in the soft X-ray band and hard X-ray detection of CTA 1. 
The absorption in the low-energy band is fixed as 2.8 X 10~ 21 
cm~ 2 following Tabic [2] The spectral behaviour of the pulsar is 
determ ined by the power-law fit of table 1 from ICaraveo et al.l 
(2010). The spectrum and the data points for the PWN are de- 
rived from the outer region (the data include ~ 5% contribution 
of the pulsar; see Table [5J- The solid, dotted, and dashed lines 
represent the uncertainty ranges of the spectral behaviour of the 
PSR, PWN, and HXD. 



Despite the exact origin of this large extended feature 
remains to be uncertain, the X-ray spectral analysis confirms 
its non-thermal nature. Assuming it is a synchrotron nebula, 
we discuss its emission properties in further details. For the 
synchrotron emitting electrons distributed as (7) <x 7~ p 
where 7 is the Lorentz factor of the wind particles, the resul- 
tant X-ray photon index F depend s on whether the em ission 
is in a fast or slow cooling regime (jCheng et al.|[2004h . Fhis 
is determined by cooling frequency of the emitting region 
v c = (18nem e c) / (orT 2 yn .B 3 ), where m e , <jt, T Byn and B are 
the electron mass, Thomson cross section, synchrotron life- 
time and the magnetic field strength respectively. We have 
r = (p + 2) /2 for a fast cooling scenario (i.e. v x > v c ) and 
r = (p + 1) /2 for slow cooling (i.e. v x < v c ). For a standard 
shock model, p spans a range of ~ 2-3 (cf. lCheng et al 1I2004J 
and references therein). This implies that V should span the 
ranges of ~ 2 — 2.5 and ~ 1.5 — 2 for the scenario of fast cool- 
ing and slow cooling, respectively. Comparing these ranges 
with the values inferred in the XIS spectrum (cf Tab. [2]), we 
suggest that the inner and the middle parts of this PWN are 
probably in a slow cooling regime. For the outskirt of the 
nebula, the best-fit photon index suggests the synchrotron 
cooling in this region is also slow. However, within the 90% 
confidence interval of the inferred photon index, it is on 
the margin of both regimes. For the inner and middle re- 
gions of the PWN, we speculate that v c should lie beyond 
the energy range covered by XIS. Assuming hv c ~ 10 keV 
and the synchrotron lifetime of the emitting electrons is 
comparable with the characteristic age of the pulsar (i.e. 
Tsyn ~ 1-4 x 10* yrs) , the magnetic field strength in the 
emitting region can be estimated at the order of B ~ 2 fiG, 
which is comparable with the typical field strength in the 
ISM (cf. iBeck et alj|2003l and references therein). 

Apart from the observation in the soft X-ray band, we 
have also investigated the field of CTA 1 with the hard X-ray 



data collected by HXD-PIN. The spectrum obtained in this 
hard band (i.e. 14 — 30 keV) can be described by a power- 
law with T = 3.33i 2 ;g*. The unabsorbed flux in 10 - 50 keV 
is found to be f x ~ (0.8 - 1.1) x 10~ n erg cm" 2 s" 1 . We 
note that the photon index inferred in this band cannot 
be tightly constrained; however, within the 90% confidence 
interval, the HXD spectrum can be smoothly connected with 
the spectrum inferred from the soft band although might 
have slightly spectral steepening as shown in Fig. [7] 

We should point out that the true nature of the hard 
X-ray emission is not yet conclusive. The FOV of the HXD 
is ~ 34' x 34', while the size of PWN in the soft X-ray is only 
< 10' x 10'. In view of the lack of imaging capability of HXD- 
PIN, the observed hard X-rays can possibly be contributed 
by other sources in the FOV of Suzaku. This might indicate 
an steeper spectral break than observed. On the other hand, 
the effects caused by the systematic uncertainty of the NXB 
background could be large. For example, if the NXB level 
is lower by a few percent (e.g. 2%, which is the no minal 
value of the NXB uncertainty; iFukazawa et al.l 120091 ) than 
the assumed value, the true hard X-ray would be higher than 
the current detection and the spectral break would not be 
as significant as observed. Si nce there are severa l upcoming 
missions, including NuST AR jHailev et al.ll2010r i and Astro- 
H (|Takahashi et all 120101 ) . which will be capable to image 
the sky up to ~ 80 keV for the first time, the exploration of 
PWNe, including the one in CTA 1, will leap into a new era 
in the near future. 
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